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I.  INTRODUCTION 


Pulsed  chemical  lasers  can  be  characterized  as  cold-reaction  or 
chain-reaction  devices.  For  the  case  of  pulsed  HF  chemical  lasers  excited 
species  are  created  by  the  cold  reaction  F  +  HF(v)  in  the  former  and 

by  the  chain  reaction  F  +  "*  HF(v)  +  H,  H  +  F2'4  HF(v)  +  F  in  the  latter.  (The 

terms  cold  and  hot  are  used  for  the  reactions  F  +  and  H  +  F^»  respectively, 
because  the  former  is  less  energetic  than  the  latter.) 

Numerous  computer  codes  have  been  developed  that  model  many  of  the 

complex  features  of  pulsed  chemical  lasers,  e.  g. ,  multilevel-multicomponent 

1  2  .3 

media,  rotational  nonequilibrium,  and  pressure  broadening.  The  accuracy 

of  these  models  primarily  is  dependent  upon  the  accuracy  of  the  rate  coefficients 

employed  therein. 

A  simplified  closed-form  analytical  solution  has  been  obtained  for  the 

4 

performance  of  cold-reaction  HF  pulsed  lasers.  Analytic  solutions  for  the 

5  6 

performance  of  cw  cold- reaction  HF  chemical  lasers  ’  can  also  be  applied 

to  cold-reaction  pulsed  lasers  by  consideration  of  the  limiting  case  of  infinitely 

fast  lateral  diffusion.  The  development  of  simplified  analytical  models  for  the 

7 

pulsed  chain- reaction  chemical  laser  has  received  less  attention.  Skifstad 
obtained  an  analytical  solution  with  the  use  of  a  "steady- state"  approximation 
for  the  chain  reaction.  The  emphasis  was  on  obtaining  numerical  estimates 
for  chain- re  action  laser  performance;  the  parametric  nature  of  the  solution  was 
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not  examined.  In  subsequent  studies  of  the  chain-reaction  chemical  laser 
the  emphasis  was  on  numerical  codes.  For  example  the  steady-state  approxi¬ 
mation  is  used  in  Reference  8  to  simplify  a  numerical  code. 

5  6 

The  simplified  models  of  cw  cold- reaction  chemical  lasers  ’  have 
yielded  useful  scaling  laws  that  characterize  laser  performance.  It  is  desire- 
able  to  deduce  the  equivalent  scaling  laws  for  pulsed  chain- reaction  chemical 
lasers.  This  is  done  herein  by  means  of  the  steady-state  approximation 
for  the  chain  reaction  and  a  simple  two-level  vibrational  energy  model.  A 
pulsed  chain- re  action  HF  chemical  laser  is  considered,  in  order  to  fix  ideas, 
and  mean  kinetic  rate  coefficients,  based  on  Reference  9,  are  employed. 
Corrections  for  the  effect  of  multiple  vibrational  levels  are  given  in  Appen¬ 
dix  A.  The  case  ^  1,  where  the  subscript  zero  denotes  initial 

concentrations,  is  of  primary  interest"*^'  **  and  is  treated  in  the  body  of 
the  report.  The  case  (^2^2^0>  *  *s  ^scusse<^  in  Appendix  B. 


-6- 


II.  THEORY 


A  pulsed  HF  laser  is  considered.  A  mixture  of  F^-H^  and  a  diluent  is 
subjected  to  an  electrical  or  flash  discharge  at  time  t  =  0,  which  partially 
dissociates  the  F.,.  Typical  initial  conditions  are  given  in  Eq.  (29).  A  chain 
reaction  is  initiated,  which  creates  vibr ationally  excited  HF.  An  approxi¬ 
mate  solution  for  the  latter  ("pumping")  process  is  given  in  the  following 
section.  The  solution  for  the  pumping  process  is  used  in  a  simple  two 
vibrational  level  model  to  obtain  scaling  laws  for  pulsed  HF  lasers  operating 
under  a  chain  reaction. 

A.  CHAIN  REACTION 

The  creation  of  HF  by  the  chain  reaction  can  be  expressed  as 

kc  3 

F  +  H,  y  (HF  )  +  H  (cold)  (la) 

l—i  V  C 

v=0 
k  6 

H  +  F2  -  y  <HFv)h  +  F  (hot)  (lb) 

v=0 

where  k  and  k.  denote  the  overall  rate  coefficients  for  the  cold  and  hot 
c  h 

reactions,  respectively.  Reverse  reactions  are  negligible.  Henceforth  such 
symbols  as  F,  H^,  and  F^  will  denote  both  species  and  concentration  in 
moles  per  cubic  centimeter.  Rate  coefficients  are  in  units  of  cubic  centimeter 
per  mole  seconds. 


The  notation 


3  6 

^T.c  =  I  <HFVc  '  ^T.h  =  2  ‘“Vh 

v=0  v=0 


hft e  hft  +  hft , 

T  T,c  T,h 

(2b) 

is  introduced,  where  HF,^,  denotes  the  total  HF  concentration. 

The  rate 

equations  can  be  expressed 

f  =  W  -  kcH2F 

(3a) 

ITT  =  k£H2F  -  k.F-.H 

(3b) 

4H  dHF 

-r-=  =  -k  H,F  - - j7^~ 

dt  c  2  dt 

(3c) 

dF,  dHF  , 

dt  =  -khF2H  =  -  dt  ' 

(3d) 

with  the  initial  conditions  (t  =  0) 

F  =  F0  ’  F2  F2,0  ’  H2  =  H2,0 

(3e) 

H  =  HFt>c  =  HFT>h  =  0 

(3£) 

In  general  H2>0  *  Fz  o  and  FQ  «  F^ 

The  following  particular  integrals  can  be  deduced 


F  +  H  =  F, 


(4a) 


HZ  +  HFT,c  -  H2,0 


(4b) 


F2  +  HFT,h  F2,0 


(4c) 


However  a  general  solution  of  Eq.  (3)  cannot  be  obtained  in  closed  form. 
Therefore  an  approximate  solution  is  sought.  The  temperature  dependence 
of  and  is  ignored. 

In  general  FQ  is  small  compared  to  F^  Q  and  Then  for  early 

times  a  "transient"  solution  can  be  obtained  in  the  form 


F  =  F 

2  *  2,0 


H2  H2,0 


(5a) 


H  F 

#■  =  1  -  =  (1  -  R)d  -  e 

0  *  0 


■t/t 


(5b) 


where 


R  [l  +  (kcH2>0/khF2,0^ 


-1 


(5c) 


tt  (kcH2,0  +  khF2,0^  1  R^(khF2,0) 


(5d) 


Here  t^  is  the  characteristic  time  for  the  transient  solution  and  is  generally 
small  compared  to  the  time  for  the  chain  reaction  to  go  to  completion. 


For  t/tt  »  i  Eq.  (5b)  indicates  that  H  and  F  approach  steady- state 
concentrations  given  by 


(6) 


From  Eq.  (6)  it  is  implied  that  dH/dt  =  dF/dt  =  0  and  dH^/dt  =  dF^/dt  in  Eq.  (3). 
Equation  (6)  can  be  used  to  characterize  H  and  F  concentrations  during  the 
remainder  of  the  chain  reaction.  Thus  for  q-  q  (i.e.,  a  chain  reaction 

terminated  by  the  disappear ance  of  H^)  Eqs.  (3c)  and  (6)  yield 


(7) 


where  t  =  (k  RF„ 
p  c  0 


,-l  ■ 


is  the  characteristic  time  for  the  chain  (pumping) 


reaction  to  go  to  completion.  (For  F-,  Q  <  Q  the  characteristic  pumping 
time  is  t  '  =  (F,  n/H,  „)  t  .)  It  also  follows  that 

P  c.,  U  U  p 


1  hft 


HFT,c  _  HFT,h 


F  -  F 

2,0  J2 


2  H 


2,0 


H 


2,0 


H 


2,0 


H 


=  1  -  e 


-t/tp 


2,0 


(8> 


The  present  solution,  Eqs.  (6)  to  (8),  characterizes  the  chain  reaction  if 
H2  5  F2  and 


r  ERzk 

p 


4(^)o 


«  i 


(9) 


The  latter  is  satisfied  when  (F/F^Jq  is  small,  as  is  generally  the  case. 
The  present  solution  is  exact  when  F^  q  =  q  and  (F/H)^  =  k^/k^ 
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B. 


PULSED  HF  LASER  MODEL 


A  two-vibrational-level  model  is  assumed.  Upper  and  lower  lasing 
levels  are  denoted  by  HFu  and  HF^,  respectively.  It  is  assumed  that  the 
pumping  reaction  creates  HF^  only,  and  that  ^  s  q. 

From  Eq.  8  HF^  is  created  by  the  pumping  reaction  at  the  rate 


/ dHF  \  dHF„  2H,  n  -t/t 
\dt  /  —  dt  "  t 

\  /  p  p 


Collisional  deactivation  by  species  M^  can  be  expressed 


HF  +  M.  -  HF.  +  M. 
u  1  1  l 


(10) 


(11a) 


The  overall  collisional  deactivation  rate  is  then 


=  -  HF 

u 


(lib) 


(lie) 


where  t  ,  — 
cd 


is  the  characteristic  collisional  deactivation  time. 


Radiative  deactivation  is  characterized  by 
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where  g  is  the  gain 


3- 


For  small  K  the  zero-power  gain  equals 


=  2  -  r  -  2e" T  +  j-  [(4  -  t  )(2  -  r )  _  8e'T]  +  0(K2) 


(16b) 


The  zero-power  gain  is  maximum  at  time 


m 


=  rricln(rM 


K  ?!  1 


(16c) 


J. 

2 


K  =  1 


Values  of  peak  zero-power  gain  are  found  by  substitution  into  Eq.  (16a). 
For  K  <  1 


G  =  0.  307  +  0(K)  K  «  1 

m  ' 

=  0.  426  K  =  1 


(I6d) 


Lasing  is  initiated  when  the  zero-power  gain  reaches  the  threshold 
value.  For  a  Fabry-Perot  resonator  with  mirror  reflectivities  and  R^ 
and  a  mirror  separation  L,  the  threshold  gain  is  gc  =  -ln(RjR2)/2L  or 

-ln(R.R?) 

G  - - - -  (17) 

C  4ctH2  qLK(K  +  1) 


The  lasing  medium  is  saturated  when  Gc  =  0  and  is  essentially  saturated  when 

G  /G  «  1,  which  for  small  K  corresponds  to  G  «  0.  307.  Under  the 
cm  c 

latter  conditions  lasing  is  initiated  at  time 
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1 


T. 

1 


G 

c 

1  +  K 


[l  +  0<oc>] 


(i8) 


which  is  small  compared  to  the  pulse  length.  For  an  unsaturated  laser 

is  found  by  sustituting  into  Eq.  (16a). 

During  lasing  G  =  Gc>  Let  E(t)  denote  the  temporal  variation  of  output 

energy  per  unit  volume.  The  intensity  of  the  output  from  a  unit  lasing  volume 

is  dE/dt  =  g  I.  The  latter  can  be  expressed,  from  Eq.  (15a)  with  G  =  G£,  as 
°c 


(1  +  K) 
K 


d  /  E(T)  \  _  2gtcd 

dT\eH2,0/'  c 


G  I 
c 


T 

=  (1  +K)e~*\  -1  _  G  (19) 

K(l+K)-1  C 

which  is  plotted  in  Fig.  2  for  the  case  Gc  =  0.  (The  quantities  Cj,  C^,  and 
C3  in  Figs.  2  and  3  represent  corrections  for  multiple  vibrational  levels  as 
discussed  in  Appendix  A.  )  The  intensity  decreases  exponentially.  (The 
initial  transient  is  not  deduced  herein. )  For  small  K  the  decrease  of  I  with 
t  is  nearly  linear.  Lasing  ends  when  1  =  0,  which  occurs  at 


r  = 
e 


K 


In 


— 


1  +  K 

1  +  GCK(1  +  K) 


(20a) 


(G  =  0,  K<  1)  (20b) 
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For  small  K,te  is  of  the  order  of  t  The  pulse  length  is  found  from  tg  -  h 


and  is  equal  to  t  for  a  saturated  laser. 

the  pulse,  is 

The  value  of  HF 

(HFT)e  K  / 

G‘  ) 

2H2,0  1  +KV 

1  +  K  / 

-  K  fl 

+  0(Cc)] 

'  1  +  K  L1 

The  value  of  HF  given  in  Eq.  (20c)  has  been  used  in  Eq. 

the  concentration  of  HF^  during  the  lasing  process. 

The  net  output  energy  per  unit  volume  per  pulse  is 

f  t 

E=Jt  gcIdt 

l 

fl  +  K  -Kt 
€H2,0|^  K  6 

(20c) 


(21) 


For  a  saturated  laser,  i.e.,  Gc  =  0,  Eq.  (21)  becomes 

E  _  .  ln(l  +  K) 

€H2,0  "  '  K 


(22a) 


2  1 

_  K  K_  K_ 

"  2  "  3  +  4 


•  •  -(K<1) 


(22b) 


The  quantity  E/fH^  q  is  a  measure  of  the  chemical  efficiency  of  the  laser. 
Equations  (16)  through  (22)  define  the  performance  of  pulsed  chain-reaction 
lasers. 
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III.  RESULTS  AND  DISCUSSION 


The  performance  of  a  pulsed  chemical  laser  operating  on  a  chain 
reaction  depends  on  the  parameters  K  and  G^.  For  a  saturated  medium 
Gc  -»  0,  and  the  performance  depends  only  on  K.  The  dependence  of  K  and 
G  on  initial  conditions  is  deduced  herein.  Simplified  expressions  for  laser 
performance  are  given  and  discussed. 

A.  SIMPLIFIED  SCALING  LAWS 

The  characteristic  pumping  rate  t  *  and  the  characteristic  collisional 
deactivation  rate  t  can  be  expressed,  respectively,  as 


and 


k  RF„ 
c  0 


1  +  Vi/A,.' 


(23a) 


=  kcdFK,0K'‘  +  k>-1’ 


HZ 

k  jH,  n 
cd  2,0 


_  M. 

2kcdMi 


(23b) 


where  denotes  collisional  deactivation  partners  other  than  HF  and  > 

e.  g.  ,  =  He>  F,  H,  and  F^  q  .  The  rate  coefficients  are  given  in  Table  L 

The  characteristic  value  of  HF,  in  Eq.  (23b),  was  taken  to  be  (HF^,)e  as  given 
in  Eq.  (20c).  The  characteristic  values  of  F  and  H  are  given  by  Eq.  (5b). 

The  relation  Kt  *  =  t  *  can  be  expressed  in  the  form 
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T 


k(t£k  + B)  =  *2 


where 


and 


2k' 


HF 

cd 


(F/F2)0(F2/H2)0 


1  +  (khF2,0/kcH2,0) 


B  = 


H2 

ccd 


2  k 


Ek 


HF 

cd 


Mi  -  \ 
cd 


1  + 


H2 

kcd  H2,0/ 


(24a) 


(24b) 


(24c) 


The  summation  term  in  Eq.  (24c)  is  generally  negligible.  Inversion  of  Eq.  (24a) 
yields 

2K  =j[4 /  +  (B  +  02)2]1/2  +  <t>Z  -  Bj  /(I  +  B)  (25) 

which  provides  an  explicit  relation  between  K  and  the  initial  conditions.  The 
quantity  <p  can  be  considered  the  basic  parameter,  rather  than  K,  which 
describes  pulsed  chemical  laser  performance.  Corresponding  values  of  K 
and  p  are  given  in  Table  I,  and  the  variation  of  output  energy  and  pulse  length 
with  <p  is  given  in  Fig.  3. 


\ 
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It  will  be  shown  that  in  general  B  «  1,  <p  «  l.and  K  «  1  for  pulsed 
chain-reaction  chemical  lasers.  Hence  K  can  be  expressed  in  the  form 

(from  Eq.  24a) 


K  =  4  [*  -  (it +  °(if]  ♦ <  B  <26a| 

=  *[4  ~  T$+  °(t)2]^  +  0(^  B  <  <^>  <  1  (26b) 

Examination  of  Eqs.  (23b)  and  (24a)  indicates  that  for  <p  «  B  collisional  deactiva 
tion  by  HF  is  negligible,  whereas  for  <p  B,  HF  is  the  major  collisional  de¬ 
activator.  The  region  near  =  B  is  the  region  where  collisional  deactivation 
by  HF  is  of  the  order  of  the  collisional  deactivation  by  the  other  collisional 
partners  (primarily  H^).  It  follows  from  Eqs.  (20b)  and  (22b)  that  for  Gc  =  0 


=  «5^-#[‘-(t)*  +  0(t)4]  »<B  (27a) 

=  +  °(*>]  B  <  ^  <  1  (27b) 


The  error  terms  are  largest  in  the  vicinity  of  B  =  <p.  The  following  expres¬ 
sions  agree  with  the  exact  solution  to  within  about  20%. 


KT  =  =  &T-  <!><  0.5  B  (28a) 

e  fH2,0  B 

-  <p  2B  <  <f>  «  1  (28b) 

The  upper  limit  in  Eq.  (28b)  is  chosen  so  as  to  satisfy  the  requirement 
t  /tg  «  1,  as  will  be  discussed.  Equations  (28a)  and  (28b)  provide  simplified 
expressions  for  saturated  laser  performance.  The  departure  of  these  expres¬ 
sions  from  the  exact  solution  [i.  e.  ,  Eqs.  (20a)  and  (22a)]  can  be  observed  in 
Fig.  3.  It  is  shown  in  Appendix  B  that  Eq.  (28)  applies  for  (H2/F2)q>  1  as 
well  as  for  (F2/H2)q  —  * • 

B.  TYPICAL  CONDITIONS 

Typical  initial  conditions  for  pulsed  chain- reaction  chemical  lasers 


HZ,0/F2.0/He  =  1/2/7 

(K/f2)0  £  to-2 

Pq  =  1  atm 

TQ  =  300  K  (29) 

The  value  of  HF^  at  the  end  of  the  lasing  pulse  is  given  by  Eq.  (20c).  A 
typical  gas  temperature  at  the  end  of  the  pulse  is  then  (neglecting  energy  in 
vibration)  in  degrees  Kelvin 


T 

e 


(6.  5  x  104)(HFT)e 


E<CvM)i 


2.  6  x  104K(1  +  K)'1 

1  +  (F2/H2)0  +  (3/5)(He/H2j0) 


(30) 
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where  the  heat  of  reaction  is  6.  5  X  10  cal/mole,  and  specific  heat  at 
constant  volume  is  3  and  5  cal/mole-K  for  monatomic  and  diatomic  gases, 
respectively. 

For  T  =  300  K 


and 


(31a) 


B  =  0. 0089  (31b) 

where  the  rate  coefficients  in  Table  II  have  been  used  and  the  summation 

term  in  Eq.  (24c)  is  neglected,  i.  e.  ,  is  the  major  collisional  deactivator 

other  than  HF.  It  follows  that  0  <  10  *  and  that  Eqs.  (28a)  and  (28b)  describe 

the  performance  of  practical  pulsed  chain-reaction  chemical  lasers.  For 

these  cases  T  -  T„  <  400  K  such  that  the  use  of  T  =  300  K  as  a  character- 
e  0  ~ 


istic  temperature  is  a  reasonable  first  estimate.  An  improved  estimate  for 
characteristic  temperature  can  be  deduced  from  Eqs.  (30)  and  (25).  Tem¬ 
perature  effects  become  important  as  the  amount  of  diluent  is  decreased. 


c. 


CORRELATION  WITH  EXPERIMENTS  AND 


NUMERICAL  STUDIES 

Equation  (28)  provides  simple  scaling  laws  for  pulsed  chain- reaction 
chemical  lasers,  which  apply  for  (^^2^  >  1 
(Appendix  B).  These  equations  have  been  deduced  by  means  of  a  two- 
vib'rational  level  model.  The  effect  of  multiple  vibrational  levels  is  to 
introduce  additional  coefficients,  as  discussed  in  Appendix  A.  Improved 
numerical  estimates  of  laser  performance  can  be  obtained  by  evaluating  these 
coefficients  to  agree  with  available  experimental  data  and  with  numerical 
code  calculations.  The  latter  procedure  is  discussed  herein. 

In  References  10  through  12  are  presented  experimental  pulsed  chemical 


,  as  well  as  for  (F^/H^Iq  ^  1 


laser  performance  data,  which  is  summarized  in  Table  II.  The  initiation 
level  (F/F2)q  reported  in  References  10  and  12  have  been  increased  by  a  factor 
of  four  (Table  II)  in  accordance  with  the  recommendations  of  Reference  13. 
Numerical  code  calculations  reported  in  Reference  3  are  included  in  Table  II. 
The  initiation  level  (F/F2)q  in  Reference  3  has  been  increased  by  a  factor  of 
two  (Table  II)  in  order  to  compensate  for  the  fact  that  the  overall  hot-reaction 


pumping  rate  in  Reference  3  is  twice  the  value  used  herein.  That  is,  the 
product  k^  (F/F2)q  in  Eq.  24b  has  been  kept  constant.  The  value  of  (F/F2 )q 
in  Reference  11  has  also  been  increased  by  a  factor  of  two  to  account  for  the 


effect  of  the  revised  value  of  k^  on  the  evaluation  of  initiation  strength  in 
Reference  11.  The  range  of  ip  in  Table  II  is  somewhat  limited  in  that  0.09  < 


Ip  <0.  33.  Equations  (28b)  and  (A-3)  are  applicable  and  indicate  that 


E/(i/jp  )  and  Et 

H2, 0 


should  be  constant. 


The  data  in  Table  II  indicate 
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(32a) 


^  E - y=  4.  6  [(1  ±  0.  5)]  J/(cm3-atm) 

H2,0 

Etg  =  0.  11  X  10‘6  [(1  ±  0.  5)]  J-sec/cm3  (32b) 

Note  that  the  data  of  Reference  12  includes  cases  with  (H^/F^ )g  >  1,  as  well 

as  (F^/H^^  >  *  anc*  *s  correiated  to  within  about  20%.  Substitution  of  the 

mean  values  from  Eqs.  (32)  into  Eq.  (A-5)  yields  Cj  =  2.  57  and  =  1.  32. 

1  /2 

Thus  the  abscissa  (C^/C^)  <f>  in  Fig.  3  can  be  replaced  by  ijj  to  within  an 

accuracy  of  about  1%.  Substitution  of  the  latter  values  of  and  into 
Eqs.  (A-2)  and  (A-  3)  yield  for  i/j  <  0.  004 

(33a) 

(33b) 

(33c) 

(34a) 

(34b) 

(34c) 
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JH 


E  =  510  J/cm3-atm 


2,0 


p_,  t  =  2.7  X  10"  atm-sec 

Ph2,0  e 

-  3  2  3 

Et  =1.4x10  ib  J- sec/cm 
e 


and  for  0. 02  <  ip  «i 


— —  =  4.  6  J  / cm3- atm 


H 


2,0 


*  8 

p  t  =  0.  11  X  10"  /ip  atm-sec 

H2,0  6 


-  6  3 

Et  =  0. 11  X  10“  J-sec/cm 
e 


Equations  (33)  and  (34)  define  the  i  ’rformance  of  pulsed  chain- reaction  HF 

chemical  lasers  for  both  (F^/H^Jq  >  1  and  (H^/F^Jq^  1  and  are  the  major 

results  of  the  present  study.  Equation  (33)  correspond  to  weak  initiation 

(where  the  assumption  Gc  »  1  may  be  violated).  In  the  latter  regime 

E  ~  (FF2/H2)q  and  t~  1  ~  H2  Q.  Equation  (34)  correspond  to  relatively  strong 

initiation  and  is  the  regime  of  practical  interest.  In  the  latter  regime 
- 1  1/2 

E~  t”  ~  (FF_)n  and  Et  =  constant, 
e  2  0  e 

As  previously  discussed,  Eqs.  (33)  and  (34)  represent  regimes  in  which 
collisional  deactivation  by  HF  is  negligible  and  dominant,  respectively.  The 
numerical  coefficients  are  considered  first  estimates,  and  further  correla¬ 
tions  with  numerical  codes  and  experimental  data  should  result  in 
improvements . 

D.  REGION  OF  VALIDITY 

In  the  present  solution  for  the  chain  reaction  pumping  process  it  is  as¬ 
sumed  that  the  initial  transient  time  t  is  small  compared  to  the  characteristic 
pump  time  t^.  This  assumption  was  shown  to  be  valid  for  (F/F2)Q  «  1  [Eq.  (9)]. 
The  resulting  solution  for  laser  performance  is  valid  provided  the  relation 
tt/te  «  1  is  also  satisfied.  Eqs.  (5d)  and  (20a)  indicate  that  for  Gc  =  0 

R2<VV(F/F2>0 

te  ~  ln(i  +  K) 

=  0(B)  0  «  B 

=  0(0)  B  «  0  «  i 

=  0(1)  0  =  0(1) 
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1 


Hence  in  order  for  the  present  solution  to  be  consistent  it  is  necessary  that 
<f>  «  1.  The  latter  relation  provides  an  upper  bound  on  the  validity  of  the 
present  theory  [e.  g. ,  Eq.  (28b)  and  Fig.  3]. 

E.  COMPARISON  WITH  COLD-REACTION  LASER 

Initial  concentrations  in  pulsed  HF  cold-reaction  chemical  lasers  are 
characterized  by  Fq,  q,  H g,  ...  ,  where  Fq  /H^  q  «  1.  By  virtue  of  the 
latter  inequality,  the  pumping  reaction  is 


HFr 


-t/t 


=  1  -  e 


t" 1  =  k  H,  » 
p  c  2,0 


(35a) 

(35b) 


The  concentration  of  HF^,,  during  the  lasing  process,  can  be  shown  to  be 
characterized  by  Fq  so  that  the  collisional  deactivation  rate  is 


-1  =kHF 
cd  cd  *0 


2H 


1  + 


2,0 


(36) 


Hence 


cd 


k  H,  . 
c  2,0 


kHFF 
cd  F0 


2H 


1  + 


2,0 


B 


(37) 
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Since  k  /kHF  =  0(10)  and  H,  n/Fn  »  1,  K  is  large  for  pulsed  cold -reaction 
lasers.  Consideration  of  Eqs.  (35a),  (10),  and  (14)  indicates  that  Eqs.  (14)  to 
(22)  describe  pulsed  cold-reaction  laser  performance  provided  that  2HZ>  Q  is 
replaced  by  FQ  therein.  Thus  for  these  lasers  K  »  1,  (HF^,)e  =  FQ  [see 
Eqs.  (20a)  and  (35a)],  and  the  energy  output  per  pulse  from  a  saturated  laser 
is  E/eFg  =  1/2.  In  contrast  the  pulsed  chain-reaction  laser  is  characterized 


IV.  CONCLUDING  REMARKS 

Scaling  laws  for  pulsed  chain-reaction  chemical  lasers  have  been 
derived.  Equations  (33)  and  (34)  provide  simple  expressions  for  the  case  of 
saturated  HF  chain  reaction  lasers.  The  results  are  applicable  for  (H^/F^Jq 
>  1  as  well  as  >  1.  It  is  hoped  that  these  results  will  be  useful  for 

the  correlation  of  data  from  experimental  and  numerical  studies  of  pulsed 
chain-reaction  chemical  laser  performance.  Such  correlations  should  pro¬ 
vide  improved  estimates  of  the  numerical  coefficients  in  Eqs.  (33)  and  (34). 

Mean  kinetic  rates  have  been  used,  and  the  numerical  coefficients  in 
Eqs.  (33)  and  (34)  have  been  based  on  values  corresponding  to  T  =  300  K. 

The  present  approach  is  applicable  when  the  amount  of  diluent  in  the  system 
is  sufficient  to  limit  the  temperature  rise  during  the  lasing  pulse 
[e.g.  ,  Eq.  (30)].  With  decrease  in  diluent,  it  may  be  necessary  to  modify  the 
coefficients  in  Eqs.  (33)  and  (34)  to  account  for  an  increased  mean  temperature 
during  the  lasing  pulse. 


APPENDIX  A 


MULTIPLE  VIBRATIONAL  LEVELS 


The  results  of  the  two -vibrational-level  model  can  be  generalized  to 
account  for  multiple  vibrational  levels.  The  generalization,  based  on  Refer¬ 
ence  6,  is  indicated  herein  for  power-on  (I  *  0)  conditions. 

The  definitions  K  =  tC(j/t  and  T  =  t/t£tj  and  the  definition  of  0[Eq.  (24b)] 
are  retained.  The  scaling  laws  deduced  from  the  two-level  model  are  gener¬ 
alized  by  making  the  following  substitutions: 

Two  Level  Multiple  Level 


K 

— 

(C1/C2)K 

T 

- 

(C2/C  t)r 

G 

c 

- 

(c3/c2)gc 

E 

- 

E/^Cj) 

- 

(■ Cl/Cz)i,Z 

The  quantities  Cj,  C^,  and  are  defined  in  Reference  6.  In  particular 
equals  the  number  of  photons  liberated  per  lasing  molecule  formed  for  the 
case  of  a  saturated  laser  with  no  collisional  deactivation.  The  quantity  is 
the  ratio  of  the  net  collisional  deactivation  rate  in  a  saturated  laser  to  the 
value  that  would  exist  if  all  the  lasing  species  molecules  were  in  the  first 
vibrational  level.  The  quantity  provides  for  a  correction  to  the  net 


collisional  deactivation  rate  for  the  case  of  a  lack  of  optical  saturation.  For 
a  given  lasing  species  the  quantities  C^,  C^,  and  are  functions  of 
exp(-2 JTr/T),  where  J  is  the  lower  lasing  level  rotational  energy  quantum 
number  and  Tn  is  the  characteristic  rotational  temperature  (=30.  16  K  for  HF). 

XV 

These  quantities  cam  be  considered  constants  of  order  one  in  the  present 
scaling  relations. 

_  j  HF  2 

Substitution  of  Eq.  (A-l)  into  Eq.  (28)  and  taking  t  =  2  kc(J  H2  Q  0 

1/2 

into  consideration  yields  for  (C^/C^)  0  £  0.  5B 


(A-2a) 


1  C1  1 

“2,0*e^C-2- 


HF 

cd 


c]  J1 

_  e  1  0 

2  cd 


1/2 


and  for  2B  <  (Cj/C-,)  0  «  1 


€H 


2,0 


.3/2 

172  * 


(A-2b) 


(A-2c) 


(A-3a) 
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l 

H,  nt  =  -s- 
2, 0  e  2 


1 

" HF  “ 

kcd  * 


(A- 3b) 


Et 


C.  1 
e  1  1 


e  ^"^2  kHF 
^  cd 


(A-3c) 


Estimates  for  Cj.  C2>  and  C3  are  needed.  For  the  present  two -level  model 


C1  ”  C2  C3  2 


(A-4) 


Estimates  for  a  multilevel  model  can  be  obtained  with  the  use  of  the  expres¬ 
sions  in  Reference  6.  An  alternative  approach  is  to  compare  the  present 
results  with  numerical  multilevel  code  solutions  or  with  experimental  data. 
The  latter  approach  is  used  herein.  In  the  regime  wherein  Eq.  (A- 3)  are 
valid 


C 


1 


(A- 5a) 


2e  k 
C2  =  Et 

e 

Substitution  of  mean  values  from  Eq. 


HF 


(32)  into  Eq. 


(A- 5)  yields 


C1  =  2.  57 
C2=  i.32 


(A-  5b) 


(A-6a) 

(A-6b) 
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3 

where  it  has  been  assumed  that  T  =  300  K  and  e  =  44.  3  x  10  J/mole.  (The 
latter  corresponds  to  radiation  at  a  wavelength  of  2.7  x  10"^  m. )  Equation 
(A-6)  indicates  that  the  present  two-level  model  underestimates  outputs  power 
by  about  a  factor  of  5  (as  a  result  of  the  neglect  of  multiple  levels  and  partial 
inversions).  Equation  (A-6)  should  be  viewed  as  a  first  estimate.  Improved 
correlations  can  be  deduced  by  further  comparisons  with  numerical  code 
results  and  experimental  data. 
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APPENDIX  B 


CASE  H2>0>F2>0 


It  has  been  assumed  in  the  body  of  the  report  that  F2  ^  >  H 2  The 

d2  q  >  F2  q  is  considered  herein. 

In  the  latter  case  the  amount  of  reactant  is  characterized  by  F2  q. 


rather  than  H0  n,  and  the  characteristic  times  are  t'  =  (F_/H,)n  t  and 
Z,  0  p  Z  Z  U  p 


cd 


2F  kHF 
^F2,0  cd 


+  B') 


- 1 


(B-  1) 


where  K*  s  t'  ,/t'  and  B'  =  (H0/F_)n  B.  The  resultant  solution  is  the  same 
cd  p  Z  Z  U 

as  that  in  the  body  of  the  report  with  F2  q  replacing  H2  Q  and  prime 
quantities  replacing  the  corresponding  unprimed  quantities.  Note  that 
t'  =t/t'cd>  <p'  =  (H2/F2)q  <p,  and 


G'  = - * - , -  (B-2) 

2a  F2>0  K-  (I  +  K')- 

For  K#  «1  it  is  found  that  t' d  =  t.d  and  K'  =  (H2/F2)QK.  As  a  result  the 
expressions  for  output  energy  E  and  pulse  time  t  for  a  saturated  chain- 
reaction  laser  in  the  limit  K  «  1,  given  in  the  body  of  the  report  [e.  g.  , 

Eqs.  (20b),  (22b),  (27),  (28),  (33),  (34)]  are  applicable  for  Q>  F2  Q  as 
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well  as  F2  q  >  q.  This  result  is  a  consequence  of  the  fact  that  for  a 
saturated  chain- reaction  HF  laser  in  the  limit  K«  1,  the  quantities  K,  K' , 
q»  and  F^  q  appear  only  as  the  products  q  K  and  F^  qK*,  and  the 
latter  are  equal. 
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LABORATORY  OPERATIONS 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel¬ 
opments  is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans¬ 
fer,  reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high- power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo¬ 
sensitive  materials  and  sensors,  high  precision  laser  ranging,  and  the  appli¬ 
cation  of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi¬ 
conducting,  superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor¬ 
rosion  and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia¬ 
tions  in  space  on  space  systems. 
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